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ABSTRACT 



Aims. Very Long Baseline Interferometry yields physical constraints on the compact radio emission of young stellar objects. At the 
same time, such measurements can be used for precise astrometric measurements of parallaxes and proper motions. Here, we aimed 
to make the first detections of very compact radio emission from class I protostars in the Corona Australis star-forming region. 
Methods. We have used the Long Baseline Array (LB A) to observe the protostars IRS 5 and IRS 7 in the Corona Australis star- 
forming region in three separate epochs. 

Results. We report the first firm radio detection of a class I protostar using very long baseline interferometry. We have detected the 
previously known non-thermal source IRS 5b, part of a binary system. IRS 5a and IRS 5N were undetected, as were all sources in the 
IRS 7 region. 

Conclusions. These results underline the unusual nature of IRS 5b as a genuine protostar with confirmed non-thermal radio emission. 
Also, these observations highlight the potential of the LBA as a tool to provide precision astrometric measurements of individual young 
stellar objects in southern star-forming regions that are not accessible to the Very Long Baseline Array in the northern hemisphere. 
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1. Introduction 



1.1. Non-thermal radio emission from young stellar objects 

Young stellar objects (YSOs) are known to be radio emitters in 
essentially all evolutionary stages, usually grouped into classes 
0— III. The earliest stages appear as class 0/1 protostars, followed 
by disk-dominated class II objects (classical T Tauri stars), with 
essentially diskless class III sources (also known as weak- line T 
Tauri stars) representing t he final pre-main sequence stage dLadal 
Il987t lAndreetai1ll993l) . In the past, it has not always been 
^* easy to disentangle thermal and non-thermal components to this 
. "Z- emission. While many YSOs show thermal emission from ion- 
ized material in jets or at the base of outflows, only very few are 
Jh confirmed non-thermal sources. An unambiguous sign of non- 
_ ^ thermal emission is the detection of circular polarization from 
gyrosynchrotron or linear polarization from synchrotron radia- 
tion. A direct determination of high brightness temperatures is 
an unambiguous sign of non-thermal emission as well. Strong 
variability is a hint of non-thermal emissio n, as are neg ative 
spectral indices (a < -0.1, where S v oc y"; [An dre 1996), but 
their interpretation is less straightforward (Giidel 2003- 

Non-thermal emission has been detected toward only a hand- 
ful of YSOs. At radio wavelength s, circular polariza tion toward 
T Tauri stars was first fou nd bv IWhite et all (1 1992b and linear 
polarization was found by Phillip s et alJ (1 19961) . In earlier evo- 
lutionary stages, the first detection of circularly polarized emis- 
sion to ward a class I protostar was reported by Feig elson et alJ 
( 1998). The target was CrA IRS 5, and this source is essentially 
still unique in its properties today. The low number of known 



non-thermal sources among YSOs could be due to the fact that 
magnetospheric emission is expected to be easily concealed by 
free-free a bsorption by ionized material along the line of sight 
(lAndrell987l) . 

Very long baseline interferometry (VLBI) radio observations 
yield direct determinations or limits on the brightness temper- 
ature of YSOs and thus constitute an important tool in this 
field. While an obvious target, CrA IRS 5 is too far south 
{6 - -37°) to be easily observable with the NRAO Very Long 
Baseline Array. The most prominent YSO studied with VLBI 
techniques so far is the T Tau system, particularly its southern 
component which is itself a multiple syst em and where non- 
thermal radio emi ssion is clearly detected (Phillips et al. 1993; 
ISmith et all 120031; Uohnston et all 12004 lLoinard et al.ll2007allbl) 
from one of the com ponents (T Tau Sb). However, detailed in- 
frared observations (Duchene et al 1 12005b have shown that the 
radio emitting object is most likely a very embedded classical T 
Tauri star rather than a protostellar source. A more robust case 
is th at of the proto-Herbig Ae/Be obj ect EC95 in the Serpens 
core (Pre ibischll 19991: iDzib et al.ll2010b . There, non-thermal ra- 
dio emission was clearly detected in Very Long Baseline Array 
(VLB A) observations from the two components of a tight binary 
system dDzib et al.ll2.010h . However, while it is clearly a very 
young stellar system, EC95 is of intermediate mass so classify- 
ing it as Class I is somewhat improper. In an attempt to find 
more such sources, iForbrich et all (l2007al) marginally detected 
the Ophiuchus source YLW 15 in VLBI observations, but higher 
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signal-to-noise observations would be needed to confirm this de- 
tection. 



1.2. Our target: CrA IRS 5 

IRS 5 in Cor o na A ustralis was first reported when 
iTavlor & Storevl §984) discovered the embedded Coronet 
cluster close to the Herbig Ae star R CrA in near- infrared ob- 
serva tions. The distance to this r egion is ~130 pc (ICasev et al.l 
1998; for a further discussion, see iNeuhauser & Forbrichll20 08). 
IRS 5 has been found to be a close binary bv Ichen & Gr aham 
d!993l) . with a separation of ~ 0'.'6 or -78 AU in Very Large 
Telescope imaging obtained by Nisin i et al. 1 (120051) . The latter 
authors also reported a spectral type of K5-7V for the brighter 
component, IRS 5a. Finally, a third component of the system 
(IRS 5N) w as found in radio observations dForbrich et al.l 
2006, l2007d) . It is located 8'.'5 north-east of IRS 5a/b, a 
projected linear separation of ~1 100 A U. Subsequent Spitzer 
and Submillimeter Array observations ([Peterson et al.l 1201 lb 
have confirmed the classification of the (unresolved) source 
IRS 5 as a class I protostar while providing a first classification 
of IRS 5N, also as a class I protostar. Interestingly, only 
IRS 5N and not IRS 5a/b is detected in millimeter continuum 
observations with the Submillimeter Array. 

In the ce ntimeter radio regime, the Coronet cluster was first 
observed by iBrownl (1 19871) who detected radio counterparts to 
most embedded infrared sources. A multi-epoch radio study 
of the cluster by ISuters et"aH (1 19961) found that the radio flux 
density of IRS 5 is variable by a factor of 5 on timescales 
of weeks, pointing to po s sible non-thermal emission. Sub- 
sequently, Feigelson et al. (1998) showed that among the ra- 
dio counterparts in this cluster, IRS 5 stands out because it 
shows circular polarization, a clear sign of non-thermal emis- 
sion. As such, CrA IRS 5 has remained one of only very few 
comparable protostars with clear signs of non-thermal emis- 
sion . Its non-thermal nature was confirmed by subseq uent stud- 
ies dForbrich et al.l 120061 l2007ct iMiettinen et al.ll2008l) . In the 
deep X-ray study of iForbrich & Preibischl (120071) . IRS 5 ap- 
pears as a marginally resolved source which could be resolved 
into two counterparts dForbrich et al.ll20~07bb . It has since been 
show n that IRS 5a and 5b have both X-ray an d radio counter- 
parts dHamaguchi et al.ll2.008l: IChoi et al. 2008b . In a polarimet- 
ric study of the region, IChoi et al.l (2009) have shown that the 
non-thermal emission of IRS 5 is entirely due to component 
IRS 5b, the northeastern source. Since the early classification 
of the IRS 5 system as a class I protostar stil l holds when in- 
corporating Spitzer data dPeterson et al.ll20Tlt) . it is still worth 
pointing out that, intriguingly, CrA IRS 5 is essentially the only 
known class I protostar with clearly non-thermal radio emission. 
The fact that it is a close binary may offer an explanation since 
this may help clearing the line of sight of ionized material. 



2. Observations and Data Reduction 

We observed the Coronet region with the Long Baseline Array 
(LB A) on 3 occasions - March 2010, October 2010 and April 
2012 (LBA project code v329). The initial goal of these obser- 
vations was an astrometric campaign to determine the distance to 
the Coronet region, and accordingly several non-thermal sources 
thought likely to be compact (IRS 5 and IRS 7) were targeted in 
the first two epochs. Initial positio ns for IRS 5 and IRS 7 wer e 
taken from Chandra observations dForbrich & Preibi sch 2007). 
In the third epoch, only IRS 5 was targeted, for reasons explained 
below. 



In each epoch, continuum observations at 8.4 GHz (3.6 cm) 
were made using all available antennas of the LBA. In all 3 
epochs, the Australia Telescope National Facility (ATNF) ob- 
servatories at Parkes (64-m single dish), the Australia Telescope 
Compact Array (ATCA; five 22-m antennas phased), and Mo- 
pra (22-m single dish), were available, as were the Hobart (26-m 
single dish) and Ceduna (30-m single dish) antennas operated 
by the University of Tasmania. For the final observation, the 12- 
m antenna at Warkworth in New Zealand operated by Auckland 
University of Technology was also used. In each case, the ob- 
servation duration was 10 hours. Two dual-polarization bands 
of width 16 MHz were sampled with two bit precision in the fre- 
quency ranges 8.409 - 8.425 GHz and 8.425 - 8.441 GHz in aU 
observations. For the final observation, additional bands span- 
ning 8.441 - 8.457 GHz and 8.457 - 8.473 GHz were utilized 
at the three ATNF stations only (making use of an additional 
digital backend available at these stations). The data were cor- 
related at Curtin U niversity using the DiFX software correlator 
dDeller et al.l 120 1 ll) in full polarization mode (correlating both 
the parallel-hand and cross-hand products). 

The target sources were phase referenced using a calibra- 
tor separated by 2° (J 1853-3628). Additionally, a bright fringe 
finder (1718-649 or 1921-293) was observed for 5 minutes be- 
fore and after every -4-hour block on the target field. For the 
first two epochs, a cycle time of 6 minutes was used, in which 
2 minutes was spent on the calibrator and two minutes on each 
target source; the predicted 1 cr image rms for each target field 
using natural weighting was 63 yt/Jy beam -1 . For the final epoch, 
the cycle time was reduced to 4.25 minutes, of which 3 minutes 
was spent on the target and 1 .25 minutes on the calibrator. Along 
with the extra bandwidth, this extra time on a single target field 
almost doubled the final sensitivity, with a predicted lcr rms of 
33 fiJy beam -1 . 

Data reduction was performed using the AIPSQ software 
packa ge, and imaging was performed using Difmap (Shepherd 
1997). The data were amplitude calibrated using system temper- 
ature data (logged from the station where available, otherwise 
a priori values used), followed by delay calibration using the 
phase reference source. Phase and amplitude calibration were 
then performed using a model of J 1853-3628 derived from these 
observations. The absolute amplitude calibration of the data is 
expected to be correct at the 20% level - the low accuracy be- 
ing due to the limited system temperature logging. Instrumental 
polarization was not corrected, as these observations were not in- 
tended for high-accuracy polarimetry and did not include scans 
on a suitable polarization calibrator. Polarization leakage at the 
LBA antenn as may be as severe as 20%, although it is likely 
much lower dDodsonl2008l) . 

After calibration, we searched for emission over a wide area 
for each target field using IMAGR in AIPS with 4 096x4 096 pix- 
els of size 0.75 milliarcseconds (mas) to cover an area of radius 
1'.'5. Fields were imaged at the position of IRS 5a/b, IRS 5N 
and IRS 7. For IRS 5a/b, the only field in which a detection was 
made, the data were corrected to re-center the image on the de- 
tection, split, averaged in frequency and exported from AIPS in 
uvfits format for imaging in Difmap. In Difmap, 512x512 pixel 
images with a pixel size of 0.5 mas were formed using natural 
weighting for all three epochs. 



http://www.aips.nrao.edu/index.shtml 



Article number, page 2 of [5] 



A. T. Deller, J. Forbrich, & L. Loinard: A VLBI Detection of the Class I Protostar CrA IRS 5 



(a) March 2010 



-36°57'22 , .'l - 



'(b) October 2010 



(c) April 2012 



19 h 01 m 48 a 08 48 s 07 19 h 01 m 48 s 08 48 s 07 19 h 01 m 48 a 08 48 s 07 

<5 (J2000.0) 6 (J2000.0) 6 (J2000.0) 

Fig. 1. Naturally weighted 3.6 cm LBA images of CrA IRS 5 at the three observed epochs. The synthesised beam is shown in the lower left 
corner of each image. The contours show -3, +3, 5, 7, 9, and 11 tr, where cr is the noise level in each image (65, 75, and 27 //Jy in the three 
images, respectively.) The scale bar in AU assumes a distance of 130 pc. 



3. Results 

Figure Q] shows the final IRS 5a/b images from difmap for all 
three epochs in Stokes /. There is a clear (llcr) detection in 
the final epoch with a peak flux density of 330 fiJy beam -1 at a 
position of 19 h 01 m 48507899 -36°57'22'/1743. A tentative de- 
tection (5.4cr, discussed further below) with a peak flux density 
of 350 fiiy beam -1 is made in the first epoch at a position of 
19 h 01 m 48?07778 -36°57'22'.'1150, while no detection is made 
in the second epoch (5cr upper limit of 375 //Jy beam -1 ). As 
shown in Figure [2] and discussed further below, the source we 
detect in our LBA images is much closer to IRS 5b than IRS 5a. 
Coupled with the previous work summarized in the introduction 
which identified IRS 5b as a non-thermal source, this positional 
coincidence allows us to conclude that we have detected IRS 5b, 
and that IRS 5a was not detected in any epochs. 

For the clear detection of IRS 5b in the final epoch, the an- 
gular resolution is 5.1 mas x 11.0 mas, at a position angle of 
— 8 1 ?6. A single gaussian component was fit to the source using 
the AIPS task JMFIT - the fitted deconvolved size is 9.8 x 5.0 
mas at a position angle of 85°, with an integrated intensity of 
575 //Jy. This corresponds to a physical size of 1.3 x 0.7 AU at 
the distance of IRS 5. However, this fit is likely an upper limit 
to the source size, as residual phase errors (no self-calibration 
was possible due to the low flux density) are likely to broaden 
the image somewhat. No signific ant signal is seen in S tokes V; 
the level of polarization (<20%; iFeigelson et aljfl998h seen in 
earlier observations would make such a detection unlikely given 
our noise level of 28 //Jy beam -1 in Stokes V and the substantial 
instrumental polarization leakage. 

We also searched fo r emission from IRS 5 N, using the nom- 
inal VLA position from Forb rich et al.l (120061) . In the final, most 
sensitive observation, the rms noise level in our LBA images at 
this position is 28 //Jy beam -1 , similar to the level reached at the 
position of IRS 5a/b. No potential sources are seen above 5<x 
within l'/5 of the nominal position, and so we can constrain the 
peak flux density of the IRS 5N source to be <140 //Jy beam -1 
at VLBI scales for this epoch. 

Our other target in the Coronet cluster, the IRS 7 region, also 
remained undetected in the two epochs it was targeted, with a 
5<r upper limit of -350 ^Jy. Wh ile neither ISuters et al.l ([l996) 
nor iForbrich et al .] (120061 l2007cl) had previously found unam- 
biguou s signs of non-thermal emission toward protostars in this 
region, Miettin en et al.l d2008) reported the tentative detection of 
gyrosynchrotron radiation from outflow lobes in this area and 
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Fig. 2. Highest-resolution (frequency 8.5 GHz, synthe sized beam 
0'.'7x0'.'31) VLA image from Fig. 12 in lChoi et al.l d2008). The tick- 
marks are spaced at 1" intervals, corresponding to 130 pc. The position 
of the source that we detect in the third epoch of our LBA observations 
is marked by a red star. We interpret the detection as a counterpart of 
IRS 5b since IRS 5a is more than twice as distant. This conclusion is 
further strengthened when the effects of proper motion are considered. 

IChoi et all (120091) reported the detection of a low degree of cir- 
cular polarization of IRS 7 A in one out of 1 3 epochs of VLA data 
analyzed. Our non-detection could thus be due to variability. 

4. Discussion 

The detection obtained here with the LBA of IRS 5b is the first 
firm VLBI detection of a well-established class I protostar. Even 
though both components of this binary were detected in radio 
emission in earlier observations with lower angular resolution 
(using the VLA), we only detect one component, IRS 5b. We 
note that the nearby component IRS 5N also remains undetected. 
Based on the peak flux density seen in our VLBI observations, 
we can derive a lower limit of >1.4 x 10 5 K for the brightness 
temperature of IRS 5b. This lower limit is too low to rule out 
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thermal emission by itself and thus does not provide a clear sig- 
nature of non-thermal emission. Such a signature was, however, 
provided earlier with multiple reported detections of circularly 
polarized emission from this source (see above). The flux den- 
sity at which the source is detected is compatible wit h the quies- 
cent fl ux densities found by Forbr ich et al.l (120061) and lChoi et al.l 
(2008). These observations constrain the size of the emission re- 
gion to less than 1 AU at the distance of CrA, by far the best con- 
straint provided by direct imaging to date. However, this is still 
much larger than the presumed size of this low-mass star, and 
the non-thermal emission is most likely due to coronal activity 
which extends at most out to a small number of stellar radii. 

As noted above, we consider the 5.4rx peak in the Stokes / 
image of IRS 5 from the first epoch observation to be a tentative 
detection, an interpretation which is strengthened by the proper 
motion which it implies. In addition to the two position mea- 
surements of IRS 5b from our first and third observing epochs, 
we can m ake use of Very L arge Array (VLA) observations pre- 
sented by Choi et all d2008l) to greatly extend the time baseline 
of these radio observations and check for consistency. We use 
the reference p osition for IRS 5b from the observation Tr 1 in 
IChoi et ail (120081) . which offers the longest time baseline and the 
highest angular resolution (VLA A-array observations made in 
1996 at a wavelength of 3.5 cm). With these three points, we 
are able to fit for a radio proper motion, which we find to be 
p^cosc) = 6.9±0.3 mas yr _1 , /jg = -27.3+0.1 mas yr . This fit 
is dominated by the two LBA positions, since the error bars are 
much smaller, but the fit is consistent with all 3 epochs. Fur- 
thermore, it closely matches the value of (fi a cosd, fi$) = (5.5, 
-27.0) mas y r -1 obtained for t he me an proper motion of CrA 
members by iNeuhauser et al.l (l2000h . To perform this fit, the 
parallax was fixed to 7.7 mas, representing a distance of 130 pc. 
The position s and fitted motio n of IRS 5b are shown in Figure|3] 
We note that lChoi et al.l (120081) were unable to determine a sig- 
nificant proper motion from VLA observations of IRS 5b alone, 
but the LBA observations added here extend the time baseline 
and have a resolution and positional accuracy more than an or- 
der of magnitude higher. 

The nominal positional accuracy of our significant VLBI de- 
tection in the third observing epoch (as reported by the AIPS task 
JMFIT) is 0.6 mas in right ascension and 0.3 mas in declination. 
Given that the separation to the primary calibrator is 2°, sys- 
tematic astromet ric errors for LBA observ ations are unlikely to 
exceed this level dDeller et al.ll2008l [2012b . With 4 such epochs, 
it should be possible to measure the parallax and proper motion 
of IRS 5b at the level of 0.2 mas and 0.2 mas yr~' respectively, 
sufficient to obtain the distance to the source with an accuracy of 
only a few pc. However, obtaining this result will require signif- 
icant additional observing resources, as a full observing track at 
each epoch is required to obtain even a low-significance detec- 
tion, and the effect of variability is still uncertain, since we have 
only one firm detection. Over the relevant observing timescale 
(several years) the binary orbital moti on would be neg ligible, 
since the binary separation is >100 AU dChoi et alJl2008h . 

5. Summary 

This first VLBI detection of a class I protostar underlines the pe- 
culiar nature of CrA IRS 5. Compared to previous detections, 
IRS 5b shows a relatively low radio flux density in our detec- 
tions. Future observations, if coinci ding with a more active 
state of the system (like that seen by IChoi et all 120081) would 
place considerably more stringent constraints on the nature of 
the emission. Furthermore, our detection shows the potential 
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Fig. 3. Proper motion of IRS 5 b fitted from the tw o LBA detections 
and one VLA position taken from lChoi et al.l d2008l) . The parallax has 
been fixed based on the assumed distance of 130 pc, but the proper 
motion results are insensitive to small variations in this parameter. 



for VLBI observations of these sources to determine accurate 
parallax and proper motion measurements in an absolute refer- 
ence frame. Our preliminary results already indicate that the 
proper motion of IRS 5b is similar to that of other members of 
the Coronet cluster. The currently best distance determination of 
CrA is based on a full or bit solution of t he double-lined spectro- 
scopic binary TY CrA (ICasev et al.ll 19981) . yielding a distance of 
128±11 pc. However, that source is not located directly in the 
Coronet cluster, and in the past, distances of up to 170 pc have 
been estimated for CrA. Our LBA measurements demonstrate 
that distance determinations to star-forming regions, based on 
parallax measureme nts using cont i nuum emission from young 
stellar objects ( e.g., iMenten et al.1 l2007t lLoinard et al.l l2007bt 
iDzib et alj2010l) is also possible at southern declinations that are 
not accessible with the NRAO VLBA. Given the low detection 
rate of non-thermal emission in protostars, our observations also 
show the importance of carefully selecting suitable VLBI targets 
from preceding measurements at lower angular resolution. 
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